Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN KTpX style file v2.2) 



Chandra and XMM-Newton observations of NGC 6251 



in 
o 
o 

(N 



D. A. Evans 1 *, M. J. Hardcastle 1 ' 2 , J. H. Croston 1 ' 3 , D. M. Worrall 1 , and M. Birkinshaw 1 

1 Department of Physics, University of Bristol, Royal Fort, Tyndall Avenue, Bristol BS8 1TL, UK 

2 School of Physics, Astronomy & Mathematics, University of Hertfordshire, College Lane, Hatfield ALIO 9AB, UK 

3 Service d'Astrophysique, CEA Saclay L'Orme des Merisiers, 91191 Gif-sur-Yvette, France 



Accepted 2005 February 3. Received 2005 January 28; in original form 2004 October 15 



ABSTRACT 



> 

m 
oo 

(N 
O 

in 
o 



Of 

6 

CO 

C3 



We present new X-ray observations of the nucleus, jet and extended emission of the 
nearby radio galaxy NGC 6251 using the Chandra/ ACIS-S camera, together with a reanal- 
ysis of archival Chandra/ ACIS-I and XMM-Newton/EPIC data. We find that the nuclear X- 
ray spectrum is well-fitted with an absorbed power-law, and that there is tentative, but not 
highly significant, evidence for Fe Ka emission. We argue that the observed nuclear X-ray 
emission is likely to originate in a relativistic jet, based on the double-peaked nature, and our 
synchrotron self-Compton modelling, of the radio-to-X-ray spectral energy distribution. How- 
ever, we cannot rule out a contribution from an accretion flow. We resolve X-ray jet emission 
in three distinct regions, and argue in favour of a synchrotron origin for all three; inverse- 
Compton emission models are possible but require extreme parameters. We detect thermal 
emission on both galaxy and group scales, and demonstrate that hot gas can confine the jet, 
particularly if relativistic beaming is important. We show evidence that the radio lobe has 
evacuated a cavity in the X-ray-emitting gas, and suggest that the lobe is close to the plane of 
the sky, with the jet entering the lobe close to the surface nearest to the observer. 
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1 INTRODUCTION 

NGC 6251 is a nearby (z = 0.0244) Fanaroff Riley type I (FRI) 
radio galaxy with a supermassive black hole of mass (4-8) xlO 8 
M jFerrarese & Fonlll999l) . It exhibits complex morphology on 
a range of angular scales, and shows a rich variety of structure 
across the electromagnetic spectrum. Optically, NGC 6251 hosts 
a well-defined dust disk of diameter 730 pc for (Ho = 70 km s 
Mpc -1 ) iFerrarese & Ford 1 19991) and, based on velocity widths 
of allowed and forbidden lines, has been classified as a Seyfert 

2 galaxy IShuder & Osterbrocldfl98lh . NGC 6251 is the parent 
galaxy of a spectacular radio jet, which extends 4.5 arcmin (130 
kpc) north-west (NW) of the nucleus and has an opening an gle of 
7.4° JWaggett. Warner. & Baldwirjfl97l iPerlev. Bridle. & Willisl 
I l984j). On smalle r (milliarcsecond) scales, VLBI observations (e.g., 
Ijones et alll 986) have detected radio emission from an unresolved 
core, and a parsec-scale jet extending to the NW of the nuc leus. A 
detect ion of a parsec-s cale counterjet was cla imed bv lSudou et alJ 
1 2000), but disputed bv ljones & Wehrldi20o3) . who placed a lower 
limit of 128 on the jet/counterjet brightness ratio, implying an an- 
gle to the line of sight 9 < 40° . VLA spectral-line ra dio observa- 
tions of the inner regions of NGC 625 1 l Werner 2002 ) show 2 1 -cm 
neutral hydrogen absorption, corresponding to an intrinsic column 
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of Nh= (4.5^2 4) x 10 20 atoms cm -2 . At higher energies, it has 
been claimed that NGC 6251 is associated with the CGRO/EGRKT 
gamma-ray source 3EG J1621 + 8203 <Mukheriee et alJ I2002L 
iFoschini et all2004l) . 

X-r ay emission from NGC 625 1 was first detected with 
ROSAT (Birkins haw & Worralll Il993l) . and has more recently 
been imaged in de t ail w ith Chandra and XMM-Newton. 
iBirkinshaw & Worralll 1 19931) presented ROSAT PSPC observa- 
tions of NGC 6251, detecting nuclear emission and a contri- 
bution from a hot, X-ray emitting halo. They argued that the 
kilopa rsec-scale radio jet cannot be confined by gas pressure alone. 
iMack. Kerp. & Kleinlfl997al) reanalvsed the ROSAT data, and con- 
firmed the presence of an elliptically-shaped thermal X-ray halo 
extending out to ~ 100 kpc, reporting enhanced X-ray emission 
coincident with components of the radio jet on 100-kpc scales. 
X-ray emission from the in ner jet of NGC 6251 was imaged by 
lHardcastle & Worralll < 19991) using the ROSAT HRI. lKerp & Mackl 
120031) used a Chandra/ AClS-l observation of NGC 6251 to ar- 
gue that the pressure of the extended halo is sufficient to confine 
the radio jet. Their observation also confirmed the PSPC detec- 
tion of the outer X-ray components of the jet. A subsequent XMM- 
Newton EPIC observation of the halo 1 Sambruna et al. 2004) mea- 
sured a temperature for the gas of kT =1.6 keV. The XMM-Newton 
spectrum of the outer regions of the X-ray jet ruled o ut a thermal 
model for the X-radiation, and lSambruna et all {2004) argued that 
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the emission was due to beamed inverse-Compton scattering of cos- 
mic microwave background (CMB) photons to X-ray energies. 

The physical origin of the nuclear X-ray continuum of 
NGC 6251, as with other radio galaxies, is uncertain. Two popular 
interpretations for the origin of the emission are: (1) accretion phe- 
nomena in or near a disk; and (2) nonthermal processes (inverse- 
Compton or synchrotron emission) in a subparsec-scale jet. The 
first model is often supported by the detection of high equivalent 
width Fe Ka lines. The second interpretation is supported for ra- 
dio galaxies in general by the observed correlation between the lu- 
minosities of soft, unresolved X-ray emission and 5-GHz core ra- 
dio emission (e.g. JWorrall & Birkinshawll994LlCanosa et alJl999t 
lHardcastle & Worralll 19991). The nuclear spec trum of NGC 6251, 
as measured with A SCA iTurner et al] 1 1 997t) and XMM-Newton 
jGliozzi et afll2004l) . was modelled with a power law and a high 
equivalent width fluorescent Fe Ka emission line, suggesting that 
the X-ray emission originates from an accretion fl ow. However, the 
powe r-law continuum measured with BeppoSAX iGuainazzi et alj 
2003) had a somewhat flatter photon index and higher luminos- 
ity than the earlier ASCA observation in particular, with no Fe Ka 
line detected. Moreover, t he dou b le-p eaked SED measured by, for 
examp le. iGuainazzi et alj l2003h . IhcI fl999). and lChiaberge et ail 
(2003), led those authors to suggest that the dominant X-ray- 
production mechanism was synchrotron self-Compton (SSC) emis- 
sion from a jet. 

This paper is organised as follows. Section 2 contains a de- 
scription of the data and a summary of our analysis. The results 
from our analysis of the nucleus, jet, and extended emission are 
presented in Section 3, and are discussed in Section 4. We end with 
our conclusions in Section 5. All results presented in this paper use 
a cosmology in which f2 m ,o = 0.3 and Qa,o = 0.7, and Ho = 70 km 
s" 1 Mpc" 1 . At the redshift of NGC 6251 (z = 0.0244), 1 arcsec 
corresponds to 485 pc . All spectral fits include absorption through 
our Galaxy using iVn.Gai = 5.82 x 10 20 atoms cm -2 iMurphv et alj 
1996). When distinguishing between different model fits to the 
data, we present F-statistic results for ease of comparison with 
other papers, although we note that this method may be unreliable 
in such circumstances IProtassov et al.l2002t) . We adopt thresholds 
of 95 and 99.9 per cent for marginally and highly significant im- 
provements in the fit, respectively. Errors quoted in this paper are 
90 per cent confidence for one parameter of interest (i.e., Xmin + 
2.7), unless otherwise stated. 



2 OBSERVATIONS AND ANALYSIS METHODS 
2.1 Chandra 

NGC 6251 was observed on 2003 November 1 1 (ObsID 4130) with 
the S3 and S4 chips of the Chandra/ ACIS charge-coupled devide 
(CCD) camera. The observation was made in the VFAINT mode, 
using a 128-row subarray (giving a 1 x 16 arcmin field of view), in 
order to reduce the CCD frame time to 0.54 s and lessen the effect 
of pileup. The data were reprocessed using CIAO v3.0.2 with the 
CALDB v2.26 calibration database to create a new level-2 events 
file with grades 0, 2, 3, 4, 6, afterglow events preserved, and the 0.5- 
arcsec pixel randomization removed. To check for intervals of high 
particle background, light curves were extracted for the ACIS-S3 
chip, excluding the core. Inspection of the light curves revealed the 
presence of several flares which were removed from our analysis, 
reducing the total exposure time from 46.5 ks to 43.0 ks. 

The Chandra observation of the nucleus of NGC 625 1 is piled 



up. We demonstrate this using two independent methods. First, the 
counts extracted from a source-centred circle of radius 1.23 arc- 
sec (2.5 pixels) in the 0.5-7 keV energy band gave 0.35 counts per 
frame exposure time, for which a pileup fraction of ~ 10-15 per 
cent is predicted using the PIMMS software. As a second indepen- 
dent test, we estimated the spatial extent of the pileup by produc- 
ing an image consisting of 'afterglow' events, which are produced 
from residual charges from either cosmic-ray events or from piled- 
up events. We computed the numbers of afterglow events per unit 
area within two extraction regions: (1) a source-centred circle of 
radius 0.492 arcsec (1 pixel); and (2) an annulus of inner radius 
0.492 arcsec (1 pixel) and outer radius 1.23 arcsec (2.5 pixels). The 
ratio of these numbers was found to be 55 ± 6. We performed a 
similar test with 2-5 keV events in the same regions, and found the 
ratio to be 19 ± 1, showing that the effect of pileup is reduced in 
our annular extraction region. The techniques we used to measure 
the core spectrum in the presence of this pileup are discussed in 
Section l3~2~T1 

In addition to the new ACIS-S data, we made use of an ear- 
lier 40-ks ACIS-I observation taken on 2000 September 11 (Ob- 
sID 847). These data w ere discussed bv lKerp & Mackl l200l and 
Sambr una et alj (2004). Because the X-ray core and inner jet were 
placed on the chip gaps in this observation, the data are not use- 
ful in constraining the properties of the inner regions of the source, 
but they provide additional information on the large-scale jet and 
extended emission. 

Chandra images presented in this paper were made using the 
0.5-5.0 keV counts in the filtered ACIS-S dataset. For the analysis 
of extended emission, we fitted radial surface-brightness profiles 
with models convolved with the Chandra poin t spread function 
(PSF) , and extracted physical parameters ( Birkinsha w & Worralll 
1993). In the Chandra analysis of the core, the centroids of all 
circular and annular extraction regions were determined using the 
ZHTOOLS software (A. Vikhlinin, private communication). 



2.2 XMM-Newton 

The XMM-Newton ob servation of NGC 6251 <Gliozzi et all2004 
Sambr una et alj r2004) was extracted from the archive for com- 
parison with the Chandra datasets, and in particular to study the 
emission from the group environment and the large-scale jet. The 
data were reprocessed using the latest software and calibration files 
available from the XMM-Newton project, SAS version 6.0.0. The 
total exposure time was 49.3 ks (MOS 1 ), 49.4 ks (MOS2), and 41 . 1 
ks (pn). 

To examine any periods of high background, we extracted 
high-energy light curves for the MOS1, MOS2, and pn cameras, 
in the energy range 10-12 keV (MOS) and 12-14 keV (pn), us- 
ing the entire field of view in each case. The data were filtered to 
include events with the PATTERN=0 and FLAG=#XMMEA_EM 
(MOS) and #XMMEA_EP (pn ) attributes on l y. For a more direct 
comparison with the method of Gliozzi et al. (2004), we extracted 
low-energy light curves for each camera in the energy range 0.4— 
10 keV, using only the CCD upon which the source was located 
[CCDJD=1 (MOS); CCD.ID=4 (pn)], but excluding a circle of 
radius 35 arcsec centred on the source. PATTERN < 12 (MOS) 
and < 4 (pn), together with FLAG=0 attributes were selected. In- 
spection of both the high-energy and low-energy light curves (see 
Fig.0 revealed several periods of high background, which we shall 
discuss further in Section 13.2.11 We also note from Figure □ that 
there is little difference between the location and amplitudes of the 
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flares extracted using the high- and low-energy light curve meth- 
ods. 

The XMM-Newton observation of the nucleus does not suffer 
from significant pileup. The count rates for a source-centred cir- 
cular extraction region of radius 35 arcsec were found for each of 
the MOS1, MOS2, and pn cameras, and are given in Table □ to- 
gether with the nominal count rates at which pileup is expected 
to become important . In addition, we followed the prescription of 
Mole ndT& Sembavl 120031) to produce an image of diagonal bip- 
ixels in the MOS cameras (assigned PATTERN numbers 26-29), 
which are produced almost solely by the pile-up of two single-pixel 
events. No such events were found. 

Spectral analysis on the large-scale jet and extended emission 
regions was performed using scripts based on the SAS EVSELECT 
tool to extract spectra from all three cameras. We included the vi- 
gnetting correction by using the SAS task EVIGWEIGHT, and there- 
fore used on-axis response files created using the SAS tasks RM- 
FGEN and ARFGEN (for a detailed description of this method see 
lArnaud et all2002l) . For the jet regions, we used local background 
subtraction in order to remove thermal emission from the surround- 
ing hot gas. For the extended emission, h owever, it was nec essary 
to use a double subtraction technique (e.g. Arnaud et al. 2002) both 
for spectral analysis and for the radial surface brightness profiling 
used to study the spatial distribution of the emission. The double 
subtraction technique first takes account of the high level of particle 
background (which the vignetting corrections incorrectly weights 
up at large off-axis angles) by subtracting the counts from identi- 
cal regions of a background template file, scaled to account for any 
differences in the particle background level of the source and back- 
ground dat asets. For this step , we used the background template 
files o fRead & Ponman ( 2003) appropriate for each instrument and 
filter, with the vignetting correction applied in the same way as for 
the source events lists. The second step of the process accounts for 
the residual background by subtracting a local background taken 
from an outer region free of source emission. This step is necessary 
first because the level of X-ray background events may be different 
in the source and background datasets, as a result of different lev- 
els of Galactic absorption and cosmic variance, and also because 
the scaling required to correctly remove the particle background 
introduces an arbitrary scaling factor into the other components of 
the background that must be taken into ac count. For further details 
of the technique, see lArnaud etalU 2002). and for more discussion 
of our background analysis methods for XMM-Newton data, see 
ICrostonetal]<2003l) . 

The radial surface brightness profiles (obtained separately for 
each XMM-Newton camera) were analysed by fitting point-source- 
only models and convolved point-source plus /3-models to deter- 
mine the significance of the contribution from extended emission. 
We modelled the XMM-Newton PSF using the most recent CCF 
(Current Calibration File) components: files XRT1_PSF_0007.CCF, 
XRT2_PSF_0007.CCF and XRT3_PSF.0006.CCF, described in the 
document CAL-SRN-0 167-1-1, available from the XMM-Newton 
website 1 . 

Images were produced from the XMM-New ton data using 
the techniques described in ICroston et alJ 120031) to interpolate 
over the chip gaps and remove contaminating background point 
sources (identified using the Chandra/ AC1S-1 data and an adap- 
tively smoothed image of the XMM-Newton data for larger radii). 
The image was not vignetting corrected, as the effect of incorrectly 



weighting up particle events at large radii would dominate over any 
real emission. This means that structure at larger radii is revealed 
in our image, but that its surface brightness is underestimated. The 
resulting image was smoothed using Gaussian kernels to show the 
distribution of extended emission, and also adaptively smoothed to 
show compact structure in the regions of extended emission. 

2.3 Radio 

The X-ray and radio structures were compared using a number 
of radio datasets. The best low-resolution image available was the 
map provided by the 3CRR Atlas 2 , a 327-MHz Westerb ork Synthe- 
sis Rad io Telescope (WSRT) image first published by Ma ck et alJ 
This image has a resolution (FWHM of restoring circu- 
lar Gaussian) of 55 arcsec. On smaller scales, we used a variety 
of NRAO V ery Large Arra y (VLA) datasets, including the 1.4- 
GHz map of Werner (2002), 330-MHz maps taken for other pur- 
poses by one of us (Birkinshaw et al., in preparation) and datasets 
at 1.4, 1.6 and 8.4 GHz taken from the VLA archive, including 
the data used to make the i mages presented by I Jones et alJ 1 1986) 
and I Jones & Wehrld £l994). A full list of the VLA datasets used 
and their properties is given in Table [2] Where we obtained VLA 
data from the archive we reduced them in the standard manner 
within AIPS, calibrating the flux with respect to a standard source, 
usually 3C286. Maps made by others w ere reduced i n a s imilar 
way. We combined th e 1.6-GHz data of I Jones et alJ ll986l) and 
Jones & Wehrld il994l) . which are the most sensitive datasets avail- 
able to us, subtracting from the A-configuration dataset to account 
for the effects of core variability, to obtain a combined dataset with 
sensitiv ity and resolution somewhat h igher than that of the early 
maps of Perlev, Bridle. & Willis! 11984). This dataset is used in sev- 
eral of the images shown in this paper. In addition, high-resolution 
VLBI maps were made available to us (D. Jones, private commu- 
nication). We discuss these in Section l3.2.2l 



3 RESULTS 

3.1 Overview of X-ray properties of NGC 6251 

Figures|2|and|3|show the X-ray emission from NGC 6251 detected 
with Chandra/ ACIS-S and XMM-Newton, respectively, with radio 
contours overlaid. We divide the X-rays from NGC 625 1 into emis- 
sion from the nucleus (Section 3.2), emission from the jet (Sec- 
tion 3.3), and extended thermal emission (Section 3.4). The nu- 
cleus, the brightest component, is detected with both Chandra and 
XMM-Newton. Components of the jet are also detected with both 
telescopes, although our new Chandra/ ACIS-S observation is re- 
quired to give a good image of the inner 15 arcsec of the jet. Ex- 
tended thermal X-ray emission is observed on the galaxy scale with 
Chandra/ AC1S-1 and ACIS-S (Fig.|2j, and on the group scale with 
Chandra/ ACIS-l and XMM-Newton (Fig.[3) (see Section 3.4). 

3.2 Core 

3.2.1 Spectrum 
a) Chandra 

The Chandra/ ACIS-S spectral fitting was performed over an 
energy range 0.5-7 keV. We chose the lower bound of the energy 
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range to be 0.5 keV, because of the known degradation of the ACIS 
quant um efficiency energies below ~ 0.4-0.5 keV (Marsh all et alJ 
2004). Owing to the pileup known to be present, we tried several 
methods to extract and model the core spectrum. The first approach 
was to extract the spectrum using a source-centred circular extrac- 
tion region of radius 1.23 arcsec (2.5 pixels), and to model the effect 
of the strong pileup known to be present using the method o flDavisI 
( 2001) as implemented in SHERPA. 

The second approach was to follow lGambill et alJ <2003l) and 
reduce complications due to pileup by extracting the spectrum us- 
ing an annulus surrounding the source, sampling the wings of the 
PSF and excluding the most heavily piled-up region. The inner ra- 
dius of the annulus was chosen to be 0.492 arcsec (1 pixel), which 
encloses ~60 per cent of the core flux. We chose the outer radius 
to be 1.23 arcsec (2.5 pixels), which encloses ~ 90-95 per cent 
of the core flux. The corresponding point-like ARFs (ancillary re- 
sponse files) were corrected for the energy-dependent missing flux 
using software (M. Tsujimoto, private communication) that calcu- 
lates the encircled energy fraction in an annular extraction region 
of a model PSF created using ChaRT and MARX. This theoretical 
PSF was smoothed by 0.35 arcsec to take into account the broaden- 
ing of the PSF associated with aspect and pixelization. Our spectral 
fits using this method were consistent with those in which the ARF 
was simply scaled by the fraction of the PSF missing in the annu- 
lus (i.e., when the energy dependence was not taken into account). 
Since pileup could still affect the spectrum of the annulus, we com- 
pared the results of modelling the s pectrum both with and without 
implementing the pileup method of Da yisU200ll) . 

We compare both methods here. The spectral parameters ob- 
tained from the 2.5-pixel circle using the pileup-model [T = 
1.83 ± 0.06, normalization = (1.49 ± 0.04) x 10 -3 ] method were 
inconsistent with those from the annulus method [T — 1.61±0.03, 
normalization = (0.87±0.02) x 10~ 3 ]. When the effect of the mod- 
erate pileup present in the annulus was taken into account with the 
iDavisl J2001I) method, the spectral parameters [T = 1.72 ± 0.04, 
normalization = (0.97 ± 0.03) x 10 -3 ] were consistent to within 
~ 10% of those obtained when the effect of the pileup was not 
considered, but still inconsistent with the results obtained from the 
circular extraction region. 

We decided to investigate this discrepancy further by extract- 
ing and modelling a s pectrum from the fram e transfer streak, fol- 
lowing the method of Marshall et al. 1 2005) to simply correct for 
the fraction of source events that occurred in the frame transfer 
streak. We found the best-fitting photon index, T = 1 .63 j^' 28 an d 
the normalization to be (9.31±1.90) x 10 -4 . These spectral param- 
eters are consistent with those found using the annulus method, but 
again are inconsistent with those found using the circle method. We 
therefore adopted the annulus method for our spectral analysis, and 
attribute the discrepancy between the circle and annulus methods 
to difficulties with the pileup model in handling data sets with large 
pileup fractions. We note in passing that using the annulus method 
has the advantage of being able to model spectrally galaxy-scale 
thermal emission. 

With an off-source background extraction region that matched 
the region used for radial profiling (see Section 3.4) there were 
9,433 net 0.5-7 keV counts in the source-centred annulus. We 
grouped the data to 50 counts per PHA bin. We initially attempted 
to model the spectrum with a single, absorbed power-law. The fit 
was acceptable, but the intrinsic absorption column was poorly con- 
strained (A^h ~ 8 x 10 18 atoms cm -2 with a 90 per cent confi- 
dence upper limit ~ 3 x 10 20 atoms cm -2 ). This was consistent 
with the value measured from 21 cm radio observations [A?h = 



(4.5± 2 ^) x 10 2(1 atoms cm~ 2 . IWerne3 J2002l) l. We therefore de- 
cided to fix the absorbing column at its radio-measured value and 
performed all subsequent spectral fits using this intrinsic absorp- 
tion. An absorbed power law with a photon index F — 1.72 ± 0.04 
gave an acceptable fit (x 2 = 121 for 136 de grees of freedom (dof)). 
We found a significant improvement in the fit with the addition of 
a thermal component in the form of a collisionally-ionized plasma 
model (APEC in XSPEC), with a temperature kT = 0.20±0.08 keV, 
and the abundance fixed at 0.35 solar. This model gave a good fit to 
the data: x 2 = 107 for 134 dof, with the probability of achieving a 
larger F by chance <0.1 per cent. The spectrally measured num- 
ber of 0.5-5 keV counts in the thermal component between 0.492 
and 1.23 arcsec is 99 ± 28 (±46), in agreement with the number 
of counts found with the radial profiling analysis (see Section 3.4), 
where uncertainties are quoted as la (unbracketed) and 90 per cent 
(bracketed), for one interesting parameter. This fit is shown in Fig- 
ure|4| and the best-fitting spectral parameters given in Table|3| 

We co nfirmed the presence of the thermal emission seen by 
iGliozzi et al]<2004 using XMM-Newton data by extracting a spec- 
trum using an annular extraction region of inner radius 0.492 arcsec 
and outer radius 35 arcsec (the same as used by G liozzi et all2004|) . 
excluding contaminating emission from the jet and frame transfer 
streak. There were 11,319 net counts in the source extraction re- 
gion. The best-fitting parameters of the absorbed power-law model 
are consistent with those found using the annulus of inner radius 
0.492 arcsec and outer radius 1.23 arcsec, and the temperature of 
the thermal component on thi s scale (kT = 0.59 ± 0.05 keV) is 
consistent with that quoted by G liozzi et alJ 120041) . The measured 
number of counts in the thermal component between 0.492 and 35 
arcsec is 912 ± 132 (±185), with uncertainties quoted as above. 
This number agrees, within the errors, with the number of thermal 
counts found in the radial profile analysis. 

We s earched for the broa d 6.4-keV Fe Kq emission line re- 
ported bv lGliozzi et alJ 120041) in the spectrum extracted from the 
annulus with inner radius 0.492 arcsec and outer radius 1.23 arc- 
sec. An absorbed power law was fitted to the data over an energy 
range of 2.0-6.0 keV and 6.7-7.0 keV, and the interpolated con- 
tinuum subtracted from the data in the energy range 6.0-6.7 keV. 
The continuum parameters were consistent with those found previ- 
ously, and provided a good fit to the data (x 2 = 8.34 for 6 dof in 
the energy range 5.5-7.0 keV). A Gaussian Fe emission line, with 
a centroid allowed to vary between 6.2 and 6.7 keV, and free line 
width was then added to the model. The line energy was poorly 
constrained, and the decrease in the fit statistic (A^ 2 = 6.18 for 
three additional parameters) not significant on an F-test (the proba- 
bility of a larger F by chance is 21 per cent). As an additional test, 
the data were regrouped to 5 counts per PHA bin and the Cash- 
statistic used to evaluate the goodness of fit. Again, the line energy 
was poorly constrained, and the decrease in fit statistic (AC = 9.14) 
small. 

Further, a Gaussian Fe Kq emission line, with best-fitting pa- 
rameters matching those found with XMM-Newton iGliozzi et alJ 
2004), was added to the Chandra continuum model. The XMM- 
Newton-measuved emission line gave a decrease in the fit statistic 
of A^ 2 = 3.14 for three additional parameters that was not 
significant (the probability of a larger F by chance is 65 per cent). 
Finally, the line para meters were changed to match those detected 
with ASCA iTurner et alJ ^997). Again, the decrease in the fit 
statistic was not significant (Ax 2 = 2.79 for three additional 
parameters). However, we note that the non-detection of Fe Ka 
emission with Chandra is unsurprising, given the relatively low 
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effective area of the ACIS camera at 6.4 keV. 

b) XMM-Newton 

A goal of our XMM-Newton analysis was to make an inde- 
pendent evaluati on of the fluorescent Fe Ka emission from the nu- 
cleus claimed bv lGliozzi etafl <2004 . As mentioned in Section 2, 
the XMM-Newton observation of NGC 625 1 is severely contami- 
nated by background flares. Contamination from background flares 
is a serious concern in the analysis of the low surface-brightness 
jet and extended emission regions, whereas for the high signal-to- 
noise core spectra it is less important. However, it is still important 
to assess carefully the effect of the background flares on the nuclear 
spectral analysis. 

To achieve this, we adopted several good time interval (GTI)- 
filtering criteria based upon both the high-energy and low-energy 
light curves discussed in Section l2~2l and presented in FigureQ The 
criteria we used are described below, and the effective exposure and 
number of 4-10 and 5.5-7.5 keV counts within a 35 arcsec source- 
centred circle are shown in Table|4]for each GTI filtering method. 

(i) Conservative filtering based on 10-12 keV (MOS) and 12- 
14 keV (pn) light curves. The maximum allowed count rates were 
0.25 s" 1 (MOS) and 1 s" 1 (pn) 

(ii) 3cr filtering based on 10-12 keV (MOS) and 12-14 keV (pn) 
light curves. Good time intervals were created by excluding data 
after 44.6 ks (MOS) and 42.6 ks (pn) (corresponding to the last 
flare), and calculating the time periods when the remaining data are 
within ±3<t of the mean level. 

(iii) 3a filtering based on 0.4—10 keV light curves (MOS and 
pn). As above. 

The XMM-Newton spectral fitting was performed jointly for 
the MOS1, MOS2 and pn cameras for an energy range of 0.4—10 
keV, using each of the GTI-filtering methods described above. We 
extracted source spectra using a circular extraction region of radius 
35 arcsec centred on the source. Background spectra were extracted 
using a source-centred annular extraction region of inner radius 140 
arcsec and outer radius 280 arcsec for the MOS cameras, and an 
off-source rectangle for the pn camera. Only events with FLAG=0 
and PATTERN < 12 (MOS) and < 4 (pn) were included in our 
analysis. The data were grouped to 40 counts per PHA bin. 

As a representative example, we describe the fitting to the con- 
tinuum performed on the data filtered using the 3<r method based 
on the high-energy light curves data: the spectral parameters de- 
rived using the other GTI-filtering methods are entirely consistent 
with this case. We initially attempted to model the spectrum with a 
single, absorbed power-law, with the intrinsic absorption was fixed 
at 4.5 x 10 20 atoms cm -2 . The normalizations of all spectral com- 
ponents were ke pt free for e ach camera. For completeness, and fol- 
lowing jGljozz^^i] 120041) . an Fe Ka emission line at ~ 6.4 keV 
was modelled with a Gaussian function and is discussed later. The 
fit was poor: \ 2 = 1347 for 1135 dof. However, a substantial im- 
provement in the fit was obtained by adding a thermal component 
with a temperature kT = 0.60 ± 0.07 keV, and abundance fixed 
at 0.35 solar. The new fit was acceptable (x 2 = 1197 for 1131 
dof), with the probability of achieving a larger F by chance <C 0.1 
per cent. The counts spectrum and best-fitting model, together with 
contributions to \ 2 is shown in Figure|5| The best-fitting spectral 
parameters for the XMM-Newton observation are given in Table|3| 

We also used the XMM-Newton data to search for possible flu- 
orescent line emission from iron, and considered the statistical sig- 
nificance of its detection for each GTI-filtering method described 



above. An absorbed power law was fitted to the spectrum over an 
energy range of 2.0-5.0 keV and 7.5-10.0 keV, and the interpo- 
lated continuum subtracted from the data in the energy range 5.0- 
7.5 keV. The continuum parameters were consistent with those de- 
scribed above, and provided a good fit to the data. The energy range 
was then extended to cover 5.0-7.5 keV, and the continuum param- 
eters frozen. A Gaussian Fe emission line with a centroid allowed 
to vary between 6.0 and 6.9 keV and free line width was then added 
to the model, and the significance of its addition assessed with an 
_F-test. This spectral analysis was performed for the MOS1+MOS2 
cameras jointly, the pn camera only, and finally all three jointly. A 
summary of the significances is given in Table|5] the best-fitting pa- 
rameters are shown in Table|6| and counts spectra for the pn camera 
are presented in Figure|6| 

Our results show that a continuum-only model provides an 
adequate fit to the spectra, and that any Fe Ka detection is enta- 
tive, but not highly significant, even when an unconservative GTI- 
filtering method is used. Figure|S| in particular, illustrates that dif- 
ferent methods of extracting light curves for GTI filtering purposes, 
i.e., using high-energy or low-energy light curves, has no effect on 
our conclusions. We also note that the only evidence for the exis- 
tence of the Fe Kq line is in the data from the pn camera; however, 
this is perhaps not surprising, as the effective area of the pn cam- 
era is still somewhat higher than that of the MOS1+MOS2 cameras 
combined. 

For a more direct comparison with the results of Gliozzi et alJ 

(2004), we chose to investigate further the properties of the Fe Ka 
line measured in the spectrum of the pn camera. We chose the 
3a-filtering criterion based on the low-energy (0.4-10 keV) light 
curves described in Section |2~2| which most clo sely matches the 
GTI-filtering method used by Gli ozzi et alJ 120041. We created con- 
fidence contours of the line energy and line width, as shown in Fig- 
ure|7| This test shows that the line parameters are very poorly con- 
strained. Even when the line width was frozen at 0.4 keV and the 
fitting reperformed, the confidence contours of the line energy and 
normalization (Fig.Q show that the line parameters are still highly 
uncertain. 

3.2.2 Source variability 
a) X-ray 

For the XMM-Newton data, we attempted to confirm the intra- 
observation variability claimed by Gliozzi et al. (2004) by creating 
X-ray light curves for each of the MOS1, MOS2, and pn cameras, 
using a source-centred circular extraction region of radius 35 arc- 
sec. Background subtraction was applied using a source-centred an- 
nular extraction region of inner radius 140 arcsec and outer radius 
280 arcsec for the MOS cameras, and an off-source rectangle for 
the pn camera. Only events with FLAG=0, PATTERN < 4 (MOS) 
and < 12 (pn), and an energy 0.8-10 keV were included in our 
analysis. In this instance, the good time interva ls were chosen to be 
the same as those used bv lGliozzi et alJ 120041) . 

The XMM-Newton data were binned to a variety of time scales, 
and a x 2 analysis used to test the null hypothesis. Statistically sig- 
nificant variability was only found with a long bin time. With bins 
of 3 ks, we found that the variability found in the data from the 
pn camera was significant (x 2 = 22.35 for 13 dof; null hypoth- 
esis probability = 95.0 per cent). Further, a co-added light curve 
for all three cameras was produced and binned to 3 ks. We found 
the variability to be significant (x ~ 26.07 for 12 dof; null hy- 
pothesis probability = 99.0 per cent). We do, however, caution that 
times corresponding to periods of high background are present in 



6 D. A. Evans et al. 



this light curve. While the background level is significantly lower 
than the source count rate, we cannot be certain of its impact on the 
extracted light curves. Background-subtracted results for the indi- 
vidual XMM-Newton cameras, as well as the sum of their data, are 
given in Table|7| 

For the Chandra data, we attempted to search for variabil- 
ity on 3-ks timescales by extracting an ACIS-S light curve from a 
source-centred circle of radius 1.23 arcsec, with an off-source back- 
ground extraction region that matched the choice of background for 
the radial profiling (see Section 3.4). Unlike the XMM-Newton ob- 
servations, no significant variability was found on 3-ks timescales 
(X 2 = 13.37 for 16 dof), although this might be due to the lim- 
ited sensitivity of Chandra compared with XMM-Newton. We note 
that fitting a line of gradient of 8.31 x 10~ 4 cts s~ 2 to the data is 
significantly preferred (A\ 2 = 7.77 for one additional parameter) 
over a uniform source brightness. Background- subtracted results 
for Chandra are given in Table □ Finally, we constructed a light 
curve of the hardness ratio, H — S/H + S, where S is the 0.7-1.5 
keV count rate and H the 3-7 keV count rate, binned to 3 ks. The 
variability of the hardness ratio was not significant (\ 2 = 13.62 
for 16 dof). 

Finally, in Table [8] we tabulate the long-term 2-10 keV 
power-law luminosity and photon-index history of NGC 6251. 
The 2-10 keV power law luminosity has varied by a factor of ~5 
during 1991-2003, while not all of the photon-index values are 
mutually consistent. The Chandra photon index of 1.67 ± 0.06 is 
somewhat flatter than that measured with XMM-Newton, but it is 
entirely plausible that variability has occurred. 

b) Radio 

High-resolution VLBI observations of NGC 6251 (e.g., 
Ijones et aljll986h show a compact, self-absorbed core with a 
highly inverted spectrum and a parsec-scale jet. Emission from 
the core does not greatly dominate over that from the jet in the 
VLBI observations, even at the highest available frequency of 10.7 
GHz, where it is most prominent. We used VLBI observations of 
NGC 6251 (D. Jones, private communication) to estimate the flux 
densities of the core and jet at several epochs and frequencies. The 
results are given in Table|5| This shows that variability approaching 
a factor of 2 has occurred on timescales of years. 

c) Optical 

Hubble Space Telescope (HST) images of the central re- 
gions of NGC 625 1 have shown an unresolved optical core, sur- 
rounded by a dusty disk <Ferrarese&Fordlll999r). Optical core 
flux densities were measured bv lHardcastle & Worralll j 19991) and 
IChiaberee et al.l f2003) from archival HST observations. There now 
exist HST data from several epochs that we can use to search for any 
optical and UV variability. We performed photometry on multi- 
epoch data from the WFPC2 instrument (F555W and F814W fil- 
ters) and from the FOC instrument (F3 42W and F410M filters). 
Following Hardcastle & Worrall 1 1999), the source flux was de- 
termined using a circular extraction region of radius 5 pixels (0.2 
arcsec), using background from an annular region of inner radius 
5 pixe ls and outer radius 10 pixels. As noted by IChiaberee et all 
(2003), the largest source of systematic error in determining the 
flux from the WFPC2 data is the choice of background, as a re- 
sult of the dusty disk surrounding the nucleus. Unresolved nuclear 
emission in the FOC data dominates by a large factor over extended 
emission, so the systematic errors may not be as high. The mea- 
sured count rates were converted to flux densities at the central fre- 



quency of each filter, using the SYNPHOT photometry package 
and assuming a power law with a = 0.8. 

The reddening to NGC 6251 is the sum of Galactic and possi- 
bly intrinsic components. The Galactic value was calculated using 
a standard relationship between the neutral hydrogen column den- 
sity and the B — V colour excess iBurstein & Heileslll978t) . We 
estimated the Galactic B — V colour excess to be 0.061 mag, in 
fair agreement with the value calcu lated using the dust maps of 
Schlesel, Finkbeiner. & Davis 1 1998) (0.087 mag). To convert be- 
tween the B — V colour exces s and the extinction, we used the 
extinction curve of lSchildlll977l) . As an example, this gives rise to 
a B-band extinction, Ab = 0.25 ma g. In addition, if the hydrogen 
column that obscures the radio core I Werner 2002}) also covers the 
optical core, we estimate a further extinction by assuming the same 
extinction ratio as our Galaxy. 

The measured core flux densities for the multi-epoch HST data 
are tabulated in Table 1101 assuming both Galactic and intrinsic 
absorption. This shows that the flux densities measured with the 
WFPC2 instrument have varied by < 10 per cent over a period of 
~ 1 year. However, data from the FOC instrument are separated by 
~ 5 years, and show variability by a factor ~ 3. 

3.3 Jet 

X-ray emission from the jet of NGC 6251 comes from three re- 
gions: the 'inner jet', extending from about 1.5 arcsec to 27 arc- 
sec fro m the nucleus, which corresp onds to the ROSAT HRI detec- 
tion of lHardcastle & Worralll j 1999t> : and two regions of the outer 
jet, whi ch we call regions 1 and 2, corresponding to 'knots' 1 
and 2 of Kero & Mack 1 2003), and extending respectively between 
about 200-265 arcsec and 330^-10 arcsec from the core. The new 
Chandra/ ACIS-S observation clearly detects the inner jet and re- 
gion 1 (Fig. but region 2 falls off the CCD. The Chandra/ ACIS- 
I data give a detection of all three regions (Fig. 1 of iKerp & Mackl 
2003), though the inner jet is only weakly detected because of the 
positioning of core and jet in the area of the central chip gap. Be- 
tween the end of the inner jet and the start of region 1, (i.e., be- 
tween 30 and 200 arcsec from the core) there is a 3^kr excess 
of counts over the local background in the ACIS-S data (about 50 
0.5-5 keV counts), suggesting that there is some X-ray emission 
associated with this region of the radio jet (the 'middle jet'). The 
data are not sensitive enough to allow us to associate this detection 
with a particular component of the radio emission. In the XMM- 
Newton data both compo nents of the outer jet (region 1, discussed 
bv lSambruna et al1l2004l and region 2) are detected, although the 
spatial resolution of XMM-Newton means that the inner jet is con- 
fused with the core. For the XMM-Newton analysis of the jet and 
extended emission, we used the 'conservative' GTI filtered data de- 
scribed in Section 3.2.1b. 

The relationship between the radio jet and the inner X-ray jet 
is shown in Figure|8| The inner jet, which contains 380 ± 20 0.4- 
7.0 keV Chandra/ ACIS-S counts, appears similar to the radio emis- 
sion; for example, a bright knot ('a') at 5 arcsec from the nucleus 
is present in both radio and X-ray images, while there is a dip in 
both radio and X-ray surface brightness at around 9 arcsec ('b'). 
At around 15 arcsec from the nucleus ('c'), the X-ray jet brightens 
while the surface brightness of the radio jet drops, and there are 
hints that the X-ray jet, which is transversely resolved with Chan- 
dra, may be edge-brightened at this point; thereafter the flux in X- 
rays drops off rapidly with distance along the jet (Fig.|9j. A spec- 
tral fit to the whole jet region, using an adjacent background of the 
same shape and at the same distance from the nucleus, is well fitted 
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(X 2 = 7.6 for 15 dof) with a single power-law model absorbed with 
the Galactic column density; the photon index is 2.29i[J 13 0- 17 f° r 
one interesting parameter) and the 1-keV flux density of this region 
of the jet is 10 ± 1 njy. A thermal (MEKAL) model with moderate 
(0.3 solar) abundance is a much worse fit, with x 2 = 28.5 for 15 
dof. There is no evidence for differences in the power-law index as 
a function of position in the jet; if we divide it into inner and outer 
halves containing approximately equal numbers of counts, the fit- 
ted power-law indices for the two halves are 2. 6 ±0.3 and 2. 2 ±0.3 
respectively. 

The two outer regions of the X-ray jet correspond to brighter 
features in the low-frequency radio image, and the X-ray structure 
of region 1, at least, is well matched to the radio (Fig. [2j. There is 
no sign of knotty or compact structure in either the radio or X-ray; 
the compact feature appearing at abou t 230 arcsec from the nucleus 
in Figure 5 of Jones & We hrkliT994h does not appear in our reduc- 
tion of their VLA dataset, and we suspect it to have been an imag- 
ing artefact. We extracted spectra for region 1 from the two Chan- 
dra datasets, using rectangular regions with adjacent off-source 
background, and also from the three XMM-Newton cameras, us- 
ing a similar but somewhat broader region to account for the larger 
XMM-Newton PSF. A fit to the ACIS-S dataset alone gives a good 
fit to a power-law model with Galactic absorption and photon index 
1.96^0 21 ' DUt a joint fit to all five datasets in the 0.4-7 keV energy 
range gives a good fit (x 2 = 30.4 for 26 dof) with lower, though 
consistent, power-law index (1.68 ± 0.13) and a 1-keV flux den- 
sity of 4.7 ± 0.4 njy. We tested the significance of the detection for 
region 2 by comparing source and background counts for the three 
XMM-Newton cameras. We measured 205 background-subtracted 
0.4-7.0 keV counts in the source region, and 1686 counts in the 
background (when scaled to the size of the source region), giving a 
> 3a level detection (the detection is also significant for each sep- 
arate camera). The best-fitting power-law spectrum (x 2 = 8.2 for 7 
dof) fitted to the three XMM-Newton cameras in region 2, again us- 
ing a rectangular extraction region, has a photon index of 1.7lg 
with a 1-keV flux density of 3.4lJ'5 njy. In both these cases the 
data can be fitted equally well with thermal models with moderate 
abundance, but high temperatures are required (5 keV in the case 
of region 1, even higher in region 2). 

3.4 Extended emission 

To search for any extended emission on 10-100 kpc scales, a radial 
profile was extracted using data from the Chandra/ACIS S3 chip 
in the energy range 0.5-5.0 keV. The profile was centred on the 
nucleus, and excluded one detected background point source, as 
well as jet emission, which was excluded by ignoring a pie slice 
between the position angles 284 and 302 degrees. The radial profile 
extended from 0.98 arcsec (2 pixels) out to 84 arcsec (170 pixels), 
with background taken from an annulus of inner radius 84 arcsec 
and outer radius 118 arcsec (240 pixels) (see Figure IToi . A PSF 
was modelled using ChaRT and MARX, and was convolved with 
a Gaussian of r.m.s. width 0.35 arcsec in order to provide a good 
match with the radial profile data close to the core. 

The radial profile analysis gave an unacceptable fit to single- 
component models of either a point source or a /3-model (used to 
describe gas in hydrostatic equilibrium), but a combination of the 
two gave an acceptable fit (xt = 0.79), as shown in Figure fTTlThe 
best-fitting parameters of the /3-model were (5 — 0.5951q J345 and 
9 C = 2.527^1251 (all errors are la for two interesting parame- 
ters, i.e., Xmin + 2.3). The FWHM of the gas distribution is ~ 2 
kpc, which is comparable to that seen in other radio-loud ellipti- 



cals, such as 3C 31 (1.7 kpc. lHardcastle et alj2002^ . In our annular 
spectral extraction regions of radii 0.492 - 1.23 arcsec and 0.492 - 
35 arcsec (see Section 3.2. 1), we would expect to see 94_3g (—49 ) 
and 1499 ± 315 (±494) thermal counts respectively, where uncer- 
tainties are quoted as la (unbracketed) and 90 per cent (bracketed), 
for two interesting parameters. 

The extended emissio n on la r ger scales, previously d i scussed 
by iBirkinshaw & Worralll Jl993l) . iMack. Kerp. & KleirJ ll997ij) 
and |Sambruna et al. (2004), is clearly seen in our maps of the 
XMM-Newton data (Fig.|5J. The emission is roughly radially sym- 
metric and centred on the radio source, tho ugh there is an X-ray 
deficit to the SW of the nucleus (also seen by Ma ck. Kerp. & Klein 
Il997at) that does not appear to be related to any of the radio struc- 
ture. To characterize this emission, we extracted radial profiles for 
the three XMM-Newton cameras and for the Chondral ACIS-I data. 
The profiles were obtained in the energy range of 0.3-7 keV out 
to a radius of 450 arcsec; the inner 2 arcsec of the Chandra profile 
was excluded to avoid the chip gaps. We fitted single point-source 
models and convolved point-source plus /3-models to the XMM- 
Newton profiles. For all profiles, the point-source models were not 
an adequate fit to the data, so that the inclusion of the /3-model 
component is needed. In all cases, the inclusion of a /3-model of 
suitable parameters leads to good values for the fit statistic (an im- 
provement from x 2 = 132 for 43 dof to x 2 = 47.9 for 42 dof for 
MOS1, from x 2 = 111 for 46 dof to x 2 = 46.7 for 45 dof for 
MOS2, and from X 2 = 677 for 28 dof to x 2 = 26.5 for 27 dof 
for pn). We tested the significance of the improvement separately 
for ea ch camera using Monte Carlo simulations (as in Crost onet alj 
120041) . and in each case we found a probability of < 0.1 per cent 
that the improvement in x 2 given above could occur by chance. 
Figure \l2\ shows the pn profile with a point-source model and the 
best-fitting point-source plus /3-model. For the ACIS-I dataset, we 
included a fixed small-scale /3-model based on the ACIS-S data in 
the fit, and were then able to obtain adequate fits with an additional 
large-scale /3-model component. We find the best-fitting parame- 
ters for a joint fit to the ACIS-I, MOS1, MOS2 and pn profiles to 
be (3 — 0.901q'2 5 , r c — 116l3g arcsec (errors are la for 2 inter- 
esting parameter s), which are consistent w ith the best-fitting mod- 
els found bv lBirkinshaw & WorrallH 19931) for the ROSAT data. The 
individual datasets all give best-fitting parameters that are consis- 
tent at the joint la level. The joint fit is good (x 2 = 152 for 134 
degrees of freedom). The residuals are dominated by systematic 
uncertainties in the modelling of the XMM-Newton PSF and in the 
background subtraction at large radii. 

We extracted spectra from the three XMM-Newton cameras to 
study the extended emission identified by the radial profile analysis, 
using the 'conservative' GTI-filtering described in Section 3.2.1b. 
We used an extraction annulus between 60 and 300 arcsec and used 
the double subtraction technique described in Section 2.2 to obtain 
appropriate background spectra. The radial profile analysis shows 
that there is a significant contribution from the bright central point 
source at these radii. We included this contribution in the spectral 
fits as a fixed component with the parameters of our best-fitting core 
spectrum (an apec model of kT = 0.60 keV and Z — 0.35 solar, 
and a power-law with F — 1.85), scaling the normalizations by the 
fraction of encircled energy in the extraction region (8.3 per cent). 
In addition to this fixed component, we fitted a second APEC model 
for the extended emission. We found a best-fitting temperature of 
1.71q'3 keV and abundance of 0.4 1°'^ solar (both consistent with 
the results o flSambruna et all2004l) . with x 2 of 292 for 305 degrees 
of freedom. 

For all three XMM-Newton cameras, the total counts in the 
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(3 model between radii of 60 and 300 arcsec are consistent with 
the counts in the thermal component of the best-fitting spectrum 
described above. The total bolometric luminosity for the group at- 
mosphere, determined by integrating the /3-model, is found to be 
(7± 1) x 10 41 ergs s~\ 



4 INTERPRETATION 
4.1 Core 

The nuclear continuum spectrum of NGC 6251 detected with 
Chandra and XMM-Newton is well described by an absorbed 
power-law model. This result is consistent with two distinct mod- 
els for the nuclear X-ray emission from radio galaxies. The first is 
that the non-thermal X-radiation is dominated by emission from 
an accretion flow, via inverse-Compton upscattering of a popu- 
lation of lower-energy photons to X-ray energies. Such a model 
would be supported by the detection of broadened fluorescent Fe 
Ka emission of high equivalent width. Howev er, our analysis of the 
XMM-Newton observation first presented by iGliozzi et alj 120041) 
has shown that the detection of the Fe Ka line is not highly signifi- 
cant (a continuum-only model provides an adequate fit to the spec- 
trum), leading us to disfavour such a model. The line is not detected 
with Chandra, although this is possibly related to the relatively 
poor photon statistics of the ACIS camera at 6.4 keV. Further, the 
poorly constrained parameters using the XMM-Newton data mean 
that it is impossible to constrain its origin, much less to argue that 
it originates in the inner regions of an accretion disk. 

A second po s sible interpretation of the data 

tBirkinshaw & WorraTJ 1 19931 : IChiaberee et"ail 120031) is that 
the X-ray emission arises primarily from the subparsec radio 
VLBI jet. To investigate this model further, we calculated 
the flux density ratio of the 1-keV power-law X-ray emission 
measured with Chandra an d the parsec-scale radio jet seen in 
5-GHz VLBI observations I Jo nes et alj 119861) . We performed a 
similar calculation for the inner (kiloparsec-scale) jet, using our 
Chandra) ACIS-S observation and 8-GHz VLA observations. As 
shown in in Table 1111 the slopes and flux-density ratios of the 
X-ray and radio components of the parsec- and kiloparsec-scale 
emissions are approximately similar. This similarity also occurs 
with the nearby radio galaxy Centaurus A, where the soft nuclear 
X-ray component is interpreted as having a jet origin lEvans et alj 
2004). The similar observed variability on timescales of years of 
the radio, optical, and X-ray data might also suggest a common 
origin for the nuclear emission, although we note that the data are 
non-contemporaneous, so no firm conclusions can be drawn. 

Previous studies of the X-ray and radio flux prop- 
erties of the active nuclei of low-power radio ga laxie s 
iFabbiano et alJll984t IWorrall & Birkinshaw! 1 1994t ICanosa et ail 
1 19991 Hardcastle & Worrall! ri999h have shown a correlation be- 
tween the luminosities of soft, unresolved X-ray emission and 
5-GHz core radio emission. NGC 6251 lies near the trend- 
line established for the large B2 sample of radio galaxy cores 
JCanosaetail fT999). This is consistent with the hypothesis that 
the observed X-ray emission is associated with the parsec-scale 
VLBI jet. This is not the only possib le interpretation. For example, 
Merlo ni" Heinz. & Di Mattea 120 03) find evidence for a correlation 
between the X-ray and radio luminosities, together with the black 
hole mass, of a sample of 100 AGN, and relate the X-ray emission 
to an accretion flow. However, the majority of their target objects 
are not radio galaxies, and so may not possess relativistic parsec- 
scale jets. 



In order to place further constraints on the X-ray emis- 
sion mechanism, we constructed a spectral energy distribu- 
tion of the nucleus of NGC 6251, using observations from 
VLBI (D. Jones, private communication), ISO (Birkins haw et 
al., i n preparation), HST (p resent work; Hardcas tle & Worralll 
l200Ct IChiaberee eta?] 120031) . Cha ndra (present work) XMM- 
Newton (present wo rk), BeppoSAX iGuainazzi et al. 2003), ASCA 
iTurner et alj Il997l). ROSAT iBirkinshaw & Worralll 1 1993h . and 
CGRO/EGKET iMukheriee et alj l2002h . HST data were dered- 
dened following the method described in the preceding section. 
The SED is double-peaked, and i s similar to previously-determined 
SEDs for NGC 6251 (e.g.. lChiaberge et alJ^Ool . where the au- 
thors argued that the SED is consistent with a blazar-type model, 
with the X-rays being produced by synchrotron-self Compton 
(SSC) emission. Figure [T3l shows a plausible SSC model fit to the 
radio-to-X-ray emission of NGC 625 1 , which has been more exten- 
sively disaissed_ejsewh£rej^ 
l2003tlGhisellini. Tavecchio & Chiabergel2004 . 

It is possible, however, that in this jet-dominated scenario 
there exists a 'hidden' accretion-related component of luminosity 
of the order 10 41 ergs s _1 . Such a model would permit a fluores- 
cent Fe Ka line, and indeed for NGC 625 1 we speculate that the 
(weak) detection of the line is evidence for this model. We shall 
return to this point elsewhere (Evans et al., in preparation). It is fur- 
ther possible that the relative strength of je t- and accretion -related 
components varies over time: as noted bv IGuainazzi et all l2003h . 
states of high nuclear luminosity might correspond to periods when 
emission from a relativistic jet dominates the observed X-ray flux, 
whereas during lower-luminosity periods, emission from an accre- 
tion flow might be most prominent. This is supported by the fact 
that the most significant detection of a fluorescent Fe Ka line, with 
ASCA in 1994 jTurner et alll997l:lGuainazzi et alj2003l) . occurred 
when the continuum flux was at its lowest historical level (see Ta- 
ble[8j. 

4.2 Jet 

The natural interpretation of the inner jet is as a synchrotron 
jet similar to those seen in other FRI- type radio galaxies (e.g. 
lHardcastle. Birkinshaw. & Worralj|200ll) . The steep X-ray spec- 
trum is hard to reconcile with any other emission mechanism and, 
as seen in Table ll2l the X-ray-to-radio flux ratio is similar to those 
seen in the resolved jets of other low-power radio galaxies where 
the jet emission is interpreted as synchrotron in nature. Additional 
support for this model comes from the detection of an apparent op- 
tical counterpart of the jet <WerneJ2002t) . The optical jet is detected 
between about 8 and 20 arcsec from the nucleus, and a simple syn- 
chrotron model with a broken power-law electron energy spectrum 
can be fitted to the radio, optical and X-ray data in this region with- 
out difficulty, with the low-frequency spectral index a being 0.5 
and the high-frequency value 1.1 (Fig. ll4> . This model is similar to 
those fitted to other radio, optical and X-ray jets. 

For an equipartition magnetic field strength in this region 
(around 2 nT) the electron energy-loss timescale for an X-ray- 
emitting electron is of the order of a few hundred years, while the 
light travel time to the base of this region is ~ 1.2 x 10 4 years, ig- 
noring the effects of projection and beaming (which may be consid- 
erable, as we discuss below). The synchrotron origin of the X-ray 
emission fro m the inner jet requires in situ particle acceleration in 
this source. iKero & Mackl (2003) come to a different conclusion, 
but their calculation of the loss timescale is based on a magnetic 
field strength which is inappropriate for the inner jet.] In situ parti- 
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cle acceleration in turn requires an energy source. In th e case of the 
well-studied FRI jets in 3C 3 1 it has been argued ( Hardcas tle et alJ 
2002a ; Laing & Bridle l l2004 that the energy for the particle ac- 
celeration that produces the observed X-ray jet comes from the 
bulk deceleration of the jet flow, a process that plausibly happens 
for many other FRI jet sources in which the jet-counterjet asym- 
metry decreases with distance from the nucleus (e.g-. lLaing et alJ 
1999). If we believe that the jet-counterjet asymmetry in the case 
of NGC 6251 (whose one-sided radio jet is morphologically very 
similar to those seen in other FRI sources) is due to relativistic 
beaming, then clearly the jet does not decelerate to sub-relativistic 
speeds in the inner few kiloparsecs, as other well-studied FRI jets 
in sources of similar luminosity apparently do, but instead retains 
at least moderately relativistic speeds out to hundreds of kpc. How- 
ever, based on the radio sidedness evidence alone, we cannot rule 
out a model in which some deceleration occurs in the inner 13 kpc 
(projected) in which we see the inner X-ray jet, and in which the 
deceleration provides the energy source for the required particle 
acceleration. 

The interpretation of the outer jet is less obvious. In region 
1, the synchrotron model that we fitted to the inner jet signifi- 
cantly overpredicts the X-ray flux density; more importantly, the 
best-fitting spectral index is too flat. The two-point spectral index 
between 1.4 GHz and 1 keV is 0.95, while the best-fitting X-ray 
spectral index is 0.71, which taken at face value rules out a one- 
zone synchrotron model (the low-frequency radio spectral index is 
~ 0.5). Spectral indices as steep as 1.09 are not ruled out at the 99 
per cent confidence level, so we can devise a one-zone synchrotron 
model in which the electron energy spectrum breaks immediately 
after the radio region to a spectral index ~ 1 and then continues 
to the X-ray. We do not regard this model as impossible, but it in- 
volves ad hoc assumptions about the electron energy spectrum and 
is only marginally consistent with the X-ray data. Region 1 is much 
larger than the synchrotron loss scale for X-ray emission electrons 
in an equipartition magnetic field strength, so we have no particular 
reason to expect a one-zone synchrotron model to be applicable, but 
nevertheless it is clear that region 1 is not well described by the sort 
of model that has traditionally been fitted to the inner synchrotron 
jets of FRIs, including NGC 6251. 

Ignoring thermal models (which cannot explain the close rela- 
tionship between radio and X-ray emissio n) the alternative simple 
model for region 1 is the one favoured by Sam bruna et alJ 120041) . 
in which the jet is highly relativistic and the X-ray emission is 
caused by scattering of the boosted microwave background radi- 
ation (e.g., Tavecchio et al. 2000). This model naturally explains 
the flat spectral index, which would be expected to be similar to 
the low-frequency radio index (0.5-0.6). Such models have been 
used to describe the jets in core-dominated Fanaroff-Riley type II 
(FRII) quasars, but in general in FRIs we might expect the jets to 
be too slow to exhibit significant effects, and the model also re- 
quires small values of the angle 9 of the jet to the line of sight. 
In NGC 6251, however, we have already seen that the jet does not 
appear to decelerate to sub-relativistic speeds on 10-kpc scales, so 
that it is possible that the 100-kpc-scale jet could be produced by 
beamed inverse-Compton emission. To investigate this, we calcu- 
lated the expected inverse-Compton flux density for a grid of bulk 
Lorentz factors and angles to the line of sight 9. We assumed a 
minimum electron Lorentz factor of 10 and a maximum of 4 x 10 , 
an electron power-law index of 2 (corresponding to a spectral in- 
dex a of 0.5) and Doppler boosting appropriate to a continuous jet 
(S oc T> 2+a , where V is the Doppler factor) and took account of the 
variation of the length of the jet as a function of 9. The anisotropic 



nature of the inverse-Compton emissi on was accounted for using 
the code desc ribe d by Hard castle et alJ l2002bl) . which implements 
the results of Brunetti 1 2000). The result of our calculations was a 
grid of Lorentz factor and angle, from which we were able to inter- 
polate the range of allowed angles for a given Lorentz factor, given 
the observed X-ray flux density. The requirements for beaming pa- 
rameters are quite extreme. A bulk Lorentz factor as low as ~ 2.5 
(/3 m 0.92) is possible, but only if the jet is pointing almost directly 
along the line of sight, which is clearly not the case. For more rea- 
sonable values of 9, the required bulk Lorentz factor is larger, and 
the model requires 9 < 15° for r < 16 (J3 < 0.998) (Fig, fill. 
These high bulk speeds would have to be maintained across the 
well-resolved X-ray jet, ~ 10 kpc in width. Whereas there is some 
evidence for a spine-sheath model involving relativistic speeds on 
parsec scales in BL Lac objects (e.g. Chiab ergeetall hoOO). and 
the spine speeds in such models are comparable to those required to 
produce the observed level of X-rays through the beamed inverse- 
Compton process, we emphasise that, if this were to persist to Mpc 
scales (as required here) it would not produce the X-ray jet struc- 
ture seen in NGC 625 1 - instead, we would expect an X-ray jet that 
was significantly narrower than the radio region. 

Other constraints on the angle to the line of sight and jet 
speeds come from analysis of the jet-counterjet ratio on the par- 
sec and kiloparsec scales: for example, the lower limit on the 
high-frequency parsec-scale jet-counterjet ratio (R) determined by 
IJones & Wehrlef <2002l) . R > 128, would imply 9 < 40° if all 
the jet-counterjet asymmetry were attributed to beaming. On larger 
scales, a weak kilopars ec-scale co unterjet i s detected in some im- 
ages, implying R = 40 iPerlev. Bridle. & Willis! 1981 : if 9 < 40°, 
this would imply /3 < 0.84 on those scales. But in region 1 of the 
outer jet the jet-counterjet ratio is > 200, requiring either higher 
speeds or a change in the angle to the line of sight. All of these 
arguments are based on the assumption that the jet and counter- 
jet are steady and intrinsically symmetrical, which is dangerous, 
since some of the properties of the large-scale jet are probably in- 
fluenced by the large-scale environment, as we shall argue later. 
Although the one-sided parsec- and kiloparsec-scale jets provide 
some evidence for relativistic speeds and angles to the line of sight 
< 40°, they cannot tell us convincingly whether the constraints on 
the CMB boosting model (Fig. 1 1 5 i are reasonable. Finally, it is 
worth noting that NGC 6251 is already one of the largest known 
radio galaxies, with a projected linear size close to 2 Mpc. It is 
therefore a priori unlikely that the whole source structure is close 
to the line of sight. 

Region 2 provides few additional constraints. In this region, 
the simple synchrotron model fitted to the inner jet gives a good 
prediction of the observed 1-keV flux density (within the large er- 
rors); the weak constraints on the X-ray spectrum mean that this 
model is acceptable in that respect too. However, the large size 
of the region means that only moderate amounts of beaming are 
required to provide the observed X-ray in the boosted inverse- 
Compton model CFig. 1151 . Thus, both a synchrotron and inverse- 
Compton model are possible without difficulty for this component. 
The inverse-Compton model would require a change in jet speed 
and/or orientation, since the constraints on beaming parameters are 
not consistent with those determined for region 1. But we know 
that the jet does bend, in projection, at around region 1, so this is 
not impossible. On the other hand, the weakly detected middle jet 
(between 30 and 200 arcsec from the core) does provide an impor- 
tant additional constraint. The X-ray counts in the ACIS-S data for 
this region correspond to a 1-keV flux density, assuming a similar 
spectrum to the other components, around 1.4 nJy. But the large 



10 D. A. Evans et al. 



size of the middle jet region means that it should be a relatively 
good source of beamed IC. With parameters similar to those deter- 
mined from region 1, we overpredict the observed X-ray emission 
by a factor ~ 5. The inverse-Compton model can only explain this 
if (a) the jet speeds up significantly between the middle jet and re- 
gion 1, or (b) there is a significant change in the angle to the line of 
sight between the middle jet and region 1, in the sense that the jet 
bends towards us at region 1, with all the bending taking place per- 
pendicular to the plane of the sky; for bulk Lorentz factors around 
5-7, the flux constraints require that the middle jet have 6 > 20° . 
We regard both of these ad hoc assumptions as unlikely. 

We conclude that all possible models for the X-ray emission 
for region 1 of the jet are problematic to some extent, but that a 
synchrotron model is more consistent with what is already known 
about FRI jets and with the large-scale properties of the source. 
This then leaves us with the question of why regions 1 and 2 are 
particularly good sites for X-ray synchrotron emission and the cor- 
responding in situ particle acceleration. We return to the nature of 
these regions in the next section. 

4.3 Extended emission 

The detection of large and small-scale thermal components allows 
us to investigate the long-standing question of jet confinement in 
NGC 625 1 . We can estimate the minimum pressures in the jet com- 
ponents in the standard way, choosing a conservatively low mini- 
mum electron energy corresponding to 7 m i n = 10, assuming a low- 
energy electron energy index of 2 (which is consistent with the ob- 
served low-frequency behaviour of the jet) and using the spectrum 
that fits the radio, X-ray and optical data at higher energies. We as- 
sume no significant contribution from relativistic protons, so that 
these are true minimum-energy values. The remaining uncertain- 
ties are projection and beaming effects. Beaming in general makes 
the jet appear brighter than it really is, and so causes us to over- 
estimate the minimum pressure. Projection, which goes along with 
beaming, causes us to underestimate the true size of components, 
so again we overestimate their pressure: but it also means that they 
are further out in the thermal atmosphere than we expect. 

If we assume that the source is in the plane of the sky (mea- 
sured lengths are true lengths) and that there is no beaming, the 
results are shown in Table ll3l and plotted in Fig. l 161 We see that the 
minimum pressures are close to the external thermal pressures for 
all these regions, except for the lobe, where the internal minimum 
pressure is much lower than the external pressure, as seen in other 
sources. (Note however that the outer pressure for the lobe is an 
extrapolation of the fitted j3 model well beyond the region where 
emission is seen, so that it should be treated with caution.) As the 
effects of beaming will tend to reduce the pressures in the jet com- 
ponents, these results suggest (contrary to what was found in earlier 
work) that the jet can be close to pressure balance with the external 
medium over its entire length. We have repeated the calculations 
using two beaming models: firstly, one that can reproduce the sid- 
edness in the inner kiloparsec-scale jet (9 = 40° and (3 — 0.8), and 
secondly, one that could produce the X-ray emission from region 1 
on an inverse-Compton model (6 = 12°, f3 — 0.98). For both these 
models, the change in the ratio of radio to thermal pressure is uni- 
formly a decrease (though in the second case the thermal pressures 
at large distances are uncertain). Including the possible effects of 
beaming and projection thus only strengthens the conclusion that 
the jet is underpressured, or at best in pressure balance with the 
thermal emission at minimu m pressure. The main differen ces be- 
tween our analysis and that of Birkinsh aw & WorrallH 19931) are (1) 



the additional small-scale thermal component, (2) slightly lower 
minimum pressures in this work [as a result of slightl y more con- 
servati ve assumptions than those used by Perlev. Bri dle. & Willis! 
( 1984)] and (3) slightly higher estimated thermal pressures in the 
large-scale thermal emission. We stress that our calculation as- 
sumes that the jet is electron-positron (with no heavy content) and 
is radiating at minimum energy. 

Does the extended thermal emission tell us anything about the 
structure and dynamics of the radio source? We begin by noting that 
Fig.|3|appears to show an arm of X-ray-emitting gas at the northern 
edge of the radio lobe, si milar to the rim of m aterial seen around the 
southern lobe of 3C 449 l Hardca stle. Worrall. & B irkinsha wll998t 
ICroston et all2003h . If real, this feature would indicate that the lobe 
is behaving similarly to other FRI radio galaxies and inflating a 
large-scale cavity in the X-ray gas. The FOV of XMM-Newton pre- 
vents us from seeing whether this feature extends to the far end of 
the N lobe. More interestingly, from the point of view of the overall 
interpretation, there is a sign ificant drop in surface brightness, first 
pointed out bv iKerp & Macld 120031) . evident in Fig. [3] at roughly 
the distance of region 1 of the X-ray jet, where the radio and X-ray 
jets brighten and the extended lobe radio emission begins. (This 
drop in surface brightness is not obvious in the radial profiles of 
Fig. 1121 but these specifically exclude emission in the direction of 
the jet.) If the positional agreement between the edge of the radio 
lobe and the drop in X-ray surface brightness is not coincidental, it 
implies that the radio lobe axis is reasonably close to the plane of 
the sky - if the lobe were strongly projected, it would be hard to 
see a discontinuity in the X-ray. This is not impossible to reconcile 
with the constraints on jet angle to the line of sight 9 < 40°, since 
the lobe subtends a large angle (70°) as seen from the nucleus, pro- 
vided that the jet enters the lobe towards the surface nearer to the 
observer, though it is hard to reconcile with the very small angles 
to the line of sight required in an inverse-Compton model. If the jet 
X-ray emission is synchrotron, then the fact that there is relatively 
strong X-ray and radio emission at region 1 may be a reflection of 
particle acceleration at a pressure discontinuity, or at some other 
interaction as the jet enters the lobe. The weaker X-ray emission at 
region 2 could then be an interaction between the jet and the edge 
of the lobe, since the jet bends between region 1 and region 2, and 
must bend again if it is to connect with the 'warm spot' at the end of 
the lobe. The geometry of this model of the jet and lobe is sketched 
in Figure fTTl 

Finally, we note that, by integrating the best-fitting large-scale 
/3-model to derive the group -scale lum inosity, and using the radio- 
quiet Lx/7x relation ICroston et alj r2003). we would predict a 
temperature of ~ 0.65 keV, which is significantly less than the 
measured value. This may suggest some heating of the thermally 
emitting gas by the radio source. As the radio lobes of this source 
are large, the PdV work done in expansion is expected to be con- 
siderable. We calculate, based on the external pressure at the radius 
of the lobes, and an assumed volume of 1.5 x 10 66 m 3 (modelling 
the lobe as a sphere of radius 230 kpc, based on the extent of the 
low-frequency radio emission), that the work done by the western 
lobe is ~ 3 x 10 52 J. A simple calculation of the heat capacity of 
the entire group, C = (3/2)Nk, where N is the total number of 
particles, calculated using the spectral and spatial parameters for 
the group atmosphere, shows that if the energy transferred from 
the radio lobes is distributed throughout the group atmosphere it 
would raise the temperature of the gas by ~ 0.6 keV. Therefore, if 
the eastern lobe also contributes a similar amount of energy (it is of 
similar size), then the temperature of the NGC 625 1 group gas may 
be explained as a result of radio-source heating due to expansion of 
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the lobes. This is consistent with our o bservations of other sources 
ICroston et all2003tlCroston. Hardcastle. & Birkinshawl2005h . al- 
though the uncertainties are large. 



5 CONCLUSIONS 

We have presented results from Chondral ACIS-S and ACIS-I and 
XMM-Newton/EPIC observations of the nucleus, jet, and extended 
emission of the radio galaxy NGC 6251. For the nucleus, we find 
the following. 

(i) The X-ray spectrum is well fitted by an absorbed power-law 
with thermal emission. There is tentative, but not highly significant, 
evidence for the detection of an Fe Kq line with XMM-Newton. 

(ii) The nuclear spec tral energy distribution is double-peaked, 
similar to that found by Chiabe rge et alj 120031) . and is well-fitted 
with a synchrotron self-Compton model. This lends support to 
a model in which the nuclear X-ray emission is due to inverse- 
Compton upscattering of low energy photons in a relativistic jet to 
higher energies. 

(iii) It is plausible that some of the X-ray emission originates in 
an accretion flow (rather than a jet), and that the relative dominance 
of accretion- and jet-related components varies with time. 

For the kiloparsec-scale jet and extended emission, our con- 
clusions are as follows. 

(i) There is a well-detected inner X-ray jet with a broad-band 
spectrum and X-ray properties similar to those of other X-ray jet 
sources, implying local particle acceleration in the inner 10-20 kpc. 

(ii) Inverse-Compton models for the previously detected large- 
scale jet emission ('regions 1 and 2' in our notation) require ex- 
treme parameters; coincidences in geometry are required to explain 
the non-detection of the inner part of the jet. Although a one-zone 
synchrotron model is only marginally viable for region 1, we pre- 
fer a synchrotron over an inverse-Compton explanation for both 
this region and region 2. In the synchrotron picture, regions 1 and 
2 may be privileged sites for particle acceleration as a result of in- 
teractions between the jet and the lobe. 

(iii) Kiloparsec-scale thermal emission is present in NGC 6251, 
with properties similar to that found in other radio-loud ellipticals, 
as well as the previously known group-scale thermal emission. We 
are able to characterize the properties of both the galaxy-scale and 
group-scale components and show that the pressure of the thermal 
gas can confine the jet if it is at minimum energy and contains 
no significant contribution from heavy particles, particularly if (as 
seems likely) relativistic beaming is important. 

(iv) The detailed relationship between the northern radio lobe 
and the extended X-ray emission suggests that the radio lobe has 
evacuated a cavity in the X-ray-emitting gas, as seen in other FRI 
sources. If so, the sharpness of the boundary between gas and lobe 
suggests that the lobe is close to the plane of the sky, and that the 
jet enters the lobe close to the surface nearer the observer. 

(v) The group gas is somewhat hotter than expected for its lumi- 
nosity, and the work done in expanding the radio lobes provides a 
plausible source for the additional heat. 
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Figure 1. XMM-Newton MOS1 light curves used to evaluate good-time intervals. Left: High-energy (10-12 keV) light curve for the whole field of view. Right: 
Low-energy (0.4-10 keV) light curve for CCD 1, excluding the source. Also shown are the time/count rate filtering criteria in the conservative case (1) and the 
3<t filtering case (2 and 3), as defined in Section 3.2. 




inner jet region. ] 



Figure 2. The inner jet and region 1 of the outer jet (as defined in the text) of NGC 6251 as seen in the new Chandra/ ACIS-S data. The image is smoothed 
with a 4-arcsec FWHM circular Gaussian; superposed are contours of the 1.6-GHz radio map discussed in Section l2~3l at 0.2 X (1, 2, 4 ... ) mjy beam -1 . 
The image shown here has been rotated through —24.4 degrees. 
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Figure 3. Adaptively smoothed, combined MOS1, MOS2 and pn image in the 0.5-5 keV energy band, with radio contours overlaid. The image shows the 
relation between the X-ray and radio jet and large-scale extended emission. Radio contours are from the Westerbork 327-MHz map discussed in Section l2~3l 
Contour levels are 1,2, 4. ..16 x2 X 10 — 2 Jy beam -1 . The radio/X-ray source on the W edge of the image is presumably a background AGN, though no optical 
counterpart is detected in DSS images. 
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Figure 4. Spectral fit to data from the Chandra observation in the energy range 0.5-7 keV with a power-law and APEC model, with the abundance fixed at 
0.35 solar. The spectral extraction region was a source-centred annulus of inner radius 0.492 arcsec and outer radius 1.23 arcsec. Contributions to \ 2 are also 
shown. 



"0} LO 
o 

CD 
CO 



CO 
■*-> 

c 

D 

O 

o 

CD 
N 

"5 

E 



lO 

d 



CN 

x o 

CM 
I 

I 



t 1 1 r 





n (I Oi r i ihi:,7.;j.i 



■ i i_ . . i i . 



0.5 



1 2 5 

channel energy (keV) 




Figure 5. Spectral fit to data from the XMM-Newton observation in the energy range 0.4-10 keV with a power-law and APEC model, with the abundance fixed 
at 0.35 solar. A Gaussian line was included in the fit for completeness. Contributions to \ 2 are shown. The GTI filtering was performed using the 3cr method 
defined in Section 3.2; the continuum parameters found when using the other GTI-nTtering methods are consistent with this case. 
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Figure 6. Counts spectra for the pn camera, showing the model contributions from an absorbed power law, a Gaussian 
are the three GTI filtering methods used for the XMM-Newton nuclear spectral analysis, as defined in Section 3.2. 
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Figure 7. Left: Energy and line width confidence contours of the Fe Kq line for the pn spectrum using the 3<r GTI-filtering method based on the low-energy 
light curve. Right: Energy and normalization confidence contours of the Fe Kq line, after freezing the line width at 0.4 keV and reperforming the spectral 
fitting. This highlights the poor constraints that can be placed on the Fe Ka line-parameters, even though an unconservative GTI-filtering method was used. 
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Figure 8. The inner jet of NGC 6251. The greyscale shows Chandra X-ray counts between 0.5 and 5 keV, binned in 0.246-arcsec pixels; black is 2 counts 
pixel -1 . Superposed are contours from a radio map made from the A- and B-configuration VLA data at 1.6 GHz, with a resolution of 1.36 X 1.09 arcsec, 
at (1, 2, 4 ... ) X 0.4 mJy beam -1 . The arrows show (a) a bright knot at 5 arcsec from the nucleus visible in both radio and X-ray, (b) a dip in the surface 
brightness at both wavebands, and (c) a bright region of the X-ray jet coincident with a dip in the radio surface brightness at 15 arcsec from the nucleus. 



18 D. A. Evans et al. 




5 10 15 20 

Distance (arcsec) 



Figure 9. A profile along the inner jet of 6251 in radio (dashed line) and X-ray (solid line). Data are taken from a rectangular strip 5.7 arcsec wide, with 
adjacent background subtraction in the case of the X-ray data, and binned in 0.8-arcsec bins. The vertical dot-dashed line shows the position of the dip in the 
radio surface brightness ('b' on Fig.|SJ. The vertical lines on the X-ray points indicate the Poisson errors. 
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Figure 10. Regions used in the Chandra/ ACIS-S radial profiling analysis. Concentric circles of radii 0.98-1 18 arcsec (2-240 pixels) are shown. Background 
was extracted from the last annulus (170-240 pixels). The regions excluded from the analysis were the two large outer rectangles, the two small rectan- 
gles (masking the frame transfer streak), a small circle (masking a point source), and a pie slice between position angles 284 and 302 degrees (avoiding 
contamination from jet emission). 
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Figure 11. Background-subtracted 0.5-5 keV radial surface brightness profile for the Chondral ACIS-S data. The best-fitting model is a composite of a 
point-like component and a /3-model, with the best-fitting /3-model parameters described in the text. 
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Figure 12. Radial surface brightness profiles for the pn data in the 0.3-7 keV energy band. Left: single point-source model. Right: convolved /3-model plus 
point source, with the joint best-fitting /3-model parameters given in the text. 
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Figure 13. Radio-to-X-ray spectrum of NGC 625 1 in the source frame, with model synchrotron and SSC fits to the data assuming a spherical geometry with 
emitting radius 5 X 10~ 4 pc. The solid line discussed shows the best-fitting synchrotron model, and the dot-dashed line shows the SSC emission expected. 
The electron spectrum is assumed to extend from 7 m i n ~ 20 to 7 max ~ 4.5 X 10 5 with a number spectral index of p = 2.4, breaking by 1 at 7 ~ 600. 
The model parameters are representative only and are poorly constrained due to the historical variability of the source and non-contemporaneous nature of the 
data. The model synchrotron emission shown is truncated at a frequency 10 12 Hz, as self-absorption becomes important below this frequency. 
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Figure 14. The radio to X-ray spectrum of the inner NGC 6251 jet. Plotted are data points from the radio, an infra-red upper limit from ISO (D. Tansley, 
private communication), the optical measurements (Werner 2002) and the X-ray data from Chandra (with 1-keV flux density and spectral errors). The solid 
line is the fitted synchrotron spectrum described in the text. 




Figure 15. Constraints on the bulk Lorentz factor and angle to the line of sight for region 1 (left) and 2 (right) of the outer jet. The allowed region of parameter 
space (given the lcr range on X-ray flux density) lies between the two lines. 
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Figure 16. The thermal pressure (solid lines) and minimum radio pressure (dashed lines) in the radio components of NGC 625 1 assuming no beaming and 
a jet in the plane of the sky (see Table ll3l . The solid lines show the total pressure in the thermal component, with the two lines reflecting the lcr statistical 
uncertainties on pressure at each radius. Uncertainties on temperature are not included. 
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Figure 17. A sketch of a possible geometry for the N jet and lobe of NGC 625 1 (not to scale). The observer is on the vertical axis, and the plane of the sky 
is perpendicular to the page. Thermal emission (sketched as a single contour of constant pressure; it should be assumed to fill the entire region) provides the 
pressure to confine the jet at small distances from the nucleus. The jet enters the lobe at region 1 and the change in conditions triggers particle acceleration, 
perhaps at a shock. The jet then bends at region 1 and interacts with the edge of the lobe at region 2 before continuing to feed the distant 'warm spot' at the 
end of the lobe. The discontinuity between the lobe and the thermal X-ray emission appears relatively sharp to a distant observer in this source geometry. 
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Table 1. XMM-Newton observed 0.15-10 keV count rates in a source-centred circle of radius 35 arcsec, together with nominal count rates above which pileup 
becomes important (as defined in the XMM-Newton Users' Handbook) 



Camera 


Total exposure time (ks) 


Count rate (cts s x ) 


Pile up threshold (cts s 1 ) 


MOS1 


23.7 


0.556 ± 0.005 


0.7 


MOS2 


24.7 


0.574 ± 0.005 


0.7 


pn 


8.1 


1.973 ± 0.016 


8 



Table 2. VLA observations of NGC 625 1 used in the analysis 



Proposal 


Date 


Freq. 


Bandwidth 


Duration 


Config. 


Map resn 


Ref. 


ID 




(GHz) 


(MHz) 


(h) 




(arcsec) 




VJ38 


1985 Apr 08 


1.665 


2 X 50 


5.0 


B 


3.0 X 2.8 


1 


AK153 


1986 Aug 20 


1.452 


25 


2.8 


B 


4.9 x 3.7 




AZ032 


1987 Aug 28 


333 


64 x 0.05 


7.0 


A 


4.0 x 4.0 


4 


AZ35 


1987 Dec 27, 

1988 Jan 09 


333 


64 x 0.05 


5.0 


B 


19 x 12 




VJ49 


1988 Nov 20 


1.665 


2 x 50 


11.6 


A 


1.3 x 1.1 


2 


AP0001 


1990 Jun 27 


8.427 


7 x 3.125 


3.3 


BnA 


0.77 x 0.58 




AM322 


1991 May 09 


8.415 


2 x 50 


1.2 


D 


8.7 x 7.7 




AB746 


1995 Aug 15 


1.387 


32 x 0.4 


3.0 


A 


1.7 x 1.2 


3 



References are (l) |jonesetai]|l986l) (2) |jones & Wehrld <1994l) (3) IWerneil Eopl and Werner et al. (in prep.) (4) Birkinshaw et al, in 
prep. The resolutions quoted are major axis X minor axis (FWHM) of the elliptical Gaussian restoring beam. The duration quoted is the 
time spent on source. 



Table 3. Chandra and XMM-Newton best-fitting spectral parameters 



Observation 



Inner Extraction Radius Outer Extraction Radius Component 



Parameters 



Normalization 



Chandra/ AC1S-S 



0.492 arcsec 



0.492 arcsec 



1.23 arcsec 



35 arcsec 



Power Law 
Thermal 



Power Law 



T = 1.67 ± 0.06 
kT = 0.20 ± 0.08 keV 
Abundance = 0.35 (frozen) 
T= 1.69 ±0.09 



(9.28Io;33)xl()- 
(1.04 ± 1.22)xl0" 4 

(9.71 ±0.37)xl0~ 4 



XMM-Newton 



arcsec 



35 arcsec 



Thermal 



Power Law 



Thermal 



kT = 0.59 ± 0.05 keV 
Abundance = 0.35 (frozen) 
T= 1.88 ±0.01 



kT = 0.58 ± 0.03 keV 
Abundance = 0.35 (frozen) 



(1.81 ±0.26)xl0" 



(1.20 ±0.02)x lO- 
Cl. 22 ±0.02) x 10" 



MOS1 = 
MOS2 = 

pn = (I.27±0.02)xl0- a 
MOS 1 = (4.87 ± 3.93) x 10" 5 
MOS2 = (6.41 ± 4.06) x 10" 5 

pn = (2.07±0.28)xl0~ 4 



The thermal abundance was fixed at 0.35 solar, and the intrinsic absorption (in excess of the Galactic value) fixed at 4.50 X 10 atoms 
cm -2 . Uncertainties are 90 per cent for one parameter of interest (i.e., Xmin + The XMM-Newton parameters quoted are for the 3<r 
GTI-filtering method defined in Section 3.2. The parameters derived with the other XMM-Newton GTI-filtering methods are consistent 
with those quoted here. 



Table 4. Summary of XMM-Newton Exposures and Counts After GTI-Filtering 







MOS1 






MOS2 






pn 




GTI filtering applied 




2-10 keV 


5.5-7.5 keV 




2-10 keV 


5.5-7.5 keV 




2-10 keV 


5.5-7.5 keV 




Exposure 


Source Counts 


Source Counts 


Exposure 


Source Counts 


Source Counts 


Exposure 


Source Counts 


Source Counts 


Conservative 


26.4 ks 


4,041 


325 


27.5 ks 


4,041 


330 


15.8 ks 


6,861 


868 


3(7 (from high-energy light curves) 


43.2 ks 


6,670 


553 


43.0 ks 


6,472 


601 


38.1 ks 


18,216 


2503 


3fj (from low-energy light curves) 


43.9 ks 


6,789 


544 


43.7 ks 


6,599 


627 


39.5 ks 


19,179 


2701 
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Table 5. Significance of XMM-Newton Fe-line addition 



GTI filtering applied 


MOS1+MOS2 


pn 


MOSl+MOS2+pn 


Conservative 


100% 


16.9% 


87.2% 


3<r (from high-energy light curves) 


100% 


17.0% 


36.0% 


3cr (from low-energy light curves) 


100% 


14.8% 


9.0% 



Percentages quoted are the probability of achieving a greater F by chance. 
Table 6. Key spectral parameters after fitting Fe Ka line to the XMM-Newton data 



GTI filtering applied 


MOS1+MOS2 


pn 


MOSl+MOS2+pn 


Conservative 


r = 1.77+J5' 


15 
07 


r _i Q5+017 

i - i.»o_ 12 


p_ i O7+0.13 
1 - l.»/_ 08 




E = N/A 




E = 6.4 keV (pegged) 


E = 6.51 keV (pegged) 




cr = N/A 




ct = 0.76 ± 0.59 keV 


cr = 4.02 X 10 -8 keV (± unconstrained) 




norm = 




norm = (7.61 + ^'gi) X 10 -6 


norm = 6.52 Xl0 — 7 (± unconstrained) 


3(7 (from high-energy 


r = i.84+°; 


16 
09 


r = i-8etg;Si' 


p_i S 5+0.09 
1 -l.»5_ 004 


light curves) 


E = N/A 




E=6.83+_°; 9 J keV 


E=6.90 ± 3.29 keV 




ct = N/A 




cr = 0.56j:g;|j keV 


<t = O.Sl + ^ li keV 




norm = 




norm = (4.46t^°g X 10~ 6 ) 


norm = (2.45+_2'45 x 10 ~ 6 ) 


3a (from low-energy 


r = 1. 84+^- 


16 
09 


r-i sY+o.ii ' 


I - 1.85_ Q4 


light curves) 


E = N/A 




E=6.49 ± 0.64) 


E = 6.46°_l 3 45 




cr = N/A 




(7 = 0.38 + 0.26 


(7 = 1.77+^ ^ X 10 -2 




norm = 




norm = (3.98t°'9S x 10 ~ 6 ) 


norm = (1.60t5 eo x 10 ~ 6 ) 



Uncertainties are 90 per cent for one interesting parameter (i.e., Xmm + 2.7). 

Table 7. x 2 results for fitting a constant count rate to the XMM-Newton and Chandra data sets with 3,000 second time bins. The XMM-Newton GTI filtering 
was identical to that o flGliozzi et al]<2004l) . 



Observatory 


Camera 


Energy range (keV) 


Average Count Rate (cts s 1 ) 


X 2 1 dof 


Null hypothesis probability 


XMM-Newton 


MOS1 


0.8-10 


0.42 


8.45/13 


0.82 




MOS2 




0.43 


15.29/13 


0.29 




pn 




1.17 


22.35/13 


0.05 




MOSl+MOS2+pn 




2.03 


26.07/12 


0.01 


Chandra 


ACIS-S 


0.3-10 


0.68 


13.37/16 


0.65 



Table 8. Long-term X-ray variability history of NGC 625 1 



Date Observatory 2-10 keV PL-only luminosity Photon Index Reference 



1991 March 


ROSAT 


(1.8 ±0.5) 


xlO 42 


20 +0.6 


Birkinshaw & Worrall (1993) 


1994 October 


ASCA 


(1.7 ±0.1) 


xlO 42 


7 ,1+0.16 
^• xi -0.19 


Turner et al. (1997) 


2000 September 


Chandra 


(4.3 ± 0.4) 


xlO 42 


1.76 ±0.16 


Guainazzi et al. (2003) 


2001 July 


BeppoSAX 


(7.1 ±0.6) 


xlO 42 


1 70+ 012 


Guainazzi et al. (2003) 


2002 March 


XMM-Newton 


(5.1 ±0.1) 


xlO 42 


1.88 ±0.01 


Present work 


2003 November 


Chandra 


(6.1 ±0.2) 


xlO 42 


1.67 ±0.08 


Present work 



All luminosities quoted are unabsorbed. For ROSAT, the observed spectrum has been extrapolated to the 2-10 keV energy range. The errors 
quoted are statistical, i.e., la. 
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Table 9. VLBI flux density estimates for the radio core and parsec-scale jet at several epochs and frequencies 



Frequency 


Date 


Core flux density (Jy) 


Jet flux density (Jy) 


Total flux density (Jy) 


Reference 


1.6 GHz 


1983 March 


0.29 


0.24 


0.53 


Jones etal. (1986) 




1985 April 


0.21 ±0.03 


0.12 ±0.02 


0.33 ± 0.04 


Present work 




1988 November 


0.42 ± 0.03 


0.09 ± 0.01 


0.51 ±0.04 


Present work 




2000 May 


0.31 ± 0.02 


0.10 ±0.01 


0.41 ± 0.02 


Present work 


2.3 GHz 


1980 February 


0.25 


0.40 


0.65 


Jones etal. (1986) 


5.0 GHz 


1981 December 


0.35 


0.54 


0.89 


Jones etal. (1986) 




2000 May 


0.38 ± 0.04 


0.13 ±0.01 


0.51 ±0.05 


Present work 


8.4 GHz 


2000 May 


0.38 ± 0.04 


0.13 ±0.01 


0.51 ±0.05 


Present work 


10.7 GHz 


1978 May 


0.59 


0.30 


0.89 


Cohen & Readhead (1979). Jones et al. (1986) 


15 GHz 


2000 May 


0.39 ± 0.04 


0.11 ±0.01 


0.50 ± 0.05 


Present work 



Table 10. Multi-epoch archival HST nuclear flux densities 



Instrument Filter Date Exposure time (seconds) Pivot wavelength (A) Flux density (/xjy) 



WFPC2 


F555W 


1994 February 21 


70 


5442 


171.56 ±5.82 


WFPC2 


F555W 


1995 June 28 


800 


5442 


182.18 ± 1.77 


WFPC2 


F814W 


1995 June 28 


720 


8002 


244.95 ± 2.66 


WFPC2 


F814W 


1996 September 13 


1000 


8002 


248.38 ± 2.27 


FOC 


F342W 


1991 May 12 


1196 


3403 


34.56 ± 0.41 


FOC 


F342W 


1993 September 30 


2246 


3403 


19.36 ± 0.22 


FOC 


F342W 


1996 February 19 


1196 


3403 


54.64 ± 0.45 


FOC 


F410M 


1993 September 30 


1676 


4088 


26.70 ± 0.80 


FOC 


F410M 


1996 February 19 


1315 


4088 


91.22 ± 1.34 



The quoted values have been dereddened assuming that the optical extinction to NGC 625 1 is the sum of both the measured Galactic and 
intrinsic absorptions 



Table 11. X-ray and radio flux densities and ratios for the kiloparsec-scale and parsec-scale jets 



Component 


1 keV X-ray flux density (Jy) 


Radio flux density (Jy) (frequency) 


Ratio (xl0~ 6 ) 


Inner (kiloparsec-scale) jet 


lxlO -8 


3.1 xlO" 2 (8 GHz) 


0.32 


Core power law / pc-scale jet 


6.2xl()- 7 


~ 1 (5 GHz) 


0.62 



Table 12. X-ray-to-radio flux density ratios for a selection of FRI type radio galaxies 



Source 1-keV X-ray to 8-GHz radio Reference 

flux density ratio (x 10 -6 ) 



NGC 6251 0.32 

3C66B 0.13 

3C31 0.24 

Cen A 0.07-9.7 (various knots) 

NGC 315 0.083 

3C 296 0.069 



Present Work 

Hardcastle. Birkinshaw. & Worrall (2001) 

Hardcastle et al. (2002b) 

Hardcastle et al. (2003) 

Worra ll. Birkinshaw. & Hardcastle (2003) 

Hardcastle et al. (2005) 



The X-ray and radio measurements are over the X-ray-detected regions of the jet. 
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Table 13. Minimum and thermal pressures in radio components of NGC 6251 



Region Minimum pressure Range Inner p Outer p 

(Pa) (arcsec) (Pa) (Pa) 



Inner jet (optical) 


9.7 


X 


10" 


-13 


8-20 


4 1+0-3 


X 


io- 


-12 


11.0 


-1-0.9 
-0.6 


X 


io- 


-13 


Middle jet 


1.5 


X 


io- 


-13 


30-200 


7 O+0.9 
'■ 8 -0.8 


X 


io- 


-13 


1.3 


4-0.1 
-0.1 


X 


io- 


-13 


Outer jet region 1 


1.4 


X 


id- 


-13 


200-265 


1 O+O.l 

i - 5 -0.1 


X 


io- 


-13 


7.2 


4-1.4 
-1.5 


X 


io- 


-14 


Outer jet region 2 


3.2 


X 


ler 


-14 


330-410 


4 9+ 1 - 4 


X 


io- 


-14 


2.5 


4-1.2 
-1.3 


X 


io- 


-14 


Nlobe 


3.2 


X 


io- 


-1H 


235-1000 


9 4+1-3 


X 


io- 


-14 


2.5 


4-3.9 
-2.2 


X 


io- 


-15 



Pressure for inner jet is calculated from the count density from the combination of the small- and large-scale /3-model fits, assuming 
kT = 0.59 keV. Other pressures are calculated using the large-scale component only, assuming kT = 1.7 keV. 



